Validity of Equation-of-Motion Approach to Kondo Problem in the Large- N limit 
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The Anderson impurity model for Kondo problem is investigated for arbitrary orbit-spin degen- 
eracy N of the magnetic impurity by the equation of motion method (EOM). By employing a new 
decoupling scheme, a set self-consistent equations for the one-particle Green function are derived and 
numerically solved in the large- N approximation. For the particle-hole symmetric Anderson model 
with finite Coulomb interaction U, we show that the Kondo resonance at the impurity site exists for 
all N > 2. The approach removes the pathology in the standard EOM for N = 2, and has the same 
level of applicability as non-crossing approximation. For N — 2, an exchange field splits the Kondo 
resonance into only two peaks, consist with the result from more rigorous numerical renormalization 
group (NRG) method. The temperature dependence of the Kondo resonance peak is also discussed. 
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The Kondo effect has been a subject of intensive inves- 
tigation both experimentally and theoretically for many 
years. The Anderson impurity model Q]] has been re- 
garded as one of the successful models, which correctly 
describe the coupling between conduction electrons and 
local magnetic impurity. Although many techniques 
have been developed for solving this model J3], it is de- 
sirable to have a semi-analytic and more generic method 
that can treat finite U case at finite temperatures. It is 
believed that the effect due to finite U is important in 
determining spectroscopic, thermodynamics, and trans- 
port properties at finite temperature |0,I3]- 

The equations of motion method (EOM) HHHHEii 
might be such a candidate. This approach has been em- 
ployed to derive an analytical expression for the single 
particle Green's function of the local electron at the im- 
purity site and from which finite temperature properties 
of the system can be obtained. The EOM successfully 
yields appro ximate but correct behaviors for the resis- 
tivity lllll,ll2ll , the spin susceptibility fill E3]/ and other 
transport properties IT3l . Especially, when the temper- 
ature is above Tjc, the results based on this approach 
agree well with those from perturbative calculations. 
However, the most serious weakness in the standard 
EOM is its failure to show the Kondo resonance peak at 
the Fermi energy for symmetric Anderson model with 
finite U for spin-orbital N = 2. Recently, it has been cor- 
rectly pointed out flHl that the particle-hole symmetric 
case is the singular point for the standard EOM. How- 
ever, we note that the the non-crossing approximation 
(NCA) method with large- N expansion yj, when gen- 
eralized to the finite-?/ cases 11611 , is successful in pro- 
ducing a Kondo resonance peak at the Fermi energy for 
the particle-hole symmetric case. Therefore, it would 
be legitimate to ask whether one can develop a large- 
N EOM and the Kondo resonance peak can be recov- 
ered when N = 2. Another important issue is whether 
this kind of method can describe the effect of applied 



exchange magnetic field correctly. By addressing these 
two important issues, this Letter will establish the va- 
lidity of the EOM for treating the Kondo problem in the 
large N limit. 

We start with the Anderson impurity Hamiltonian 
with arbitrary spin-orbital degeneracy N to study the 
Kondo problem. In this model, a single band for con- 
duction electrons is adopted, and the magnetic impurity 
has N/2 degenerate localized orbitals plus each orbital 
carrying a spin degeneracy of 2. The Anderson Hamilto- 
nian of a magnetic impurity with orbit-spin degeneracy 
N in metal has the following expression: 
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where we have defined k = (k,cr) for conduction band 
indices. c\ and fl are respectively the creation opera- 
tors for conduction and / electrons at the impurity site. 
The quantities e k , e Q are the band energy of conduc- 
tion electrons and energy of the local electron in the 
a-orbital (spin index included) at the impurity site, re- 
spectively. For simplicity, a constant conduction elec- 
trons density of states is assumed, i.e., p(e) = 1/2D 
when —D<ek<D, where D is band half-width. Here an 
SU (N)-type Coulomb interaction with strength U is as- 
sumed between electrons of different orbitals (or with 
opposite spins in the same orbital) at the impurity site, 
and Vka represents s-f hybridization. 

The decoupling scheme used in the original work of 
Appelbaum, Penn, and Lacroix (APL) J5[ @1 nas been 
regarded as standard in EOM. Using this "standard" 
scheme for the EOM, one can obtain a coupled set of N 
self-consistent equations for the local electron Green's 
function. In order to make these equations close among 
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themselves in this method, truncation of the perturba- 
tive terms has to be made. However, a close look at 
the truncating approximation, certain higher order of 
Green's function terms should not be simply discarded, 
and they need to be properly included since these terms 
are absent at N = 2 and become important when N is 
larger than 2. We here generalize Lacoix's N = 2 de- 
coupling scheme to large N case and write the matrix 
element of the following higher order Green's function 
in the approximate form 

</«4/ 7 4/«'i4» = (4/«'K/a4/7i4» + o(v 2 ) 

(2) 

Corrections to this decoupling scheme are of order of 
V 2 and can be neglected in the limit of V^Q and the 
validity of this corrections has been discussed by Czy- 
choll [7]. Adding the above terms to our work is nec- 
essary to make local electron Green's functions to sat- 
isfy a set of self-consistent equations and make the EOM 
more powerful as it will be demonstrated below. After 
a complicated and lengthy derivation, the final expres- 
sion for the matrix element of impurity Green's function 
G a p = </q|4»=G„ q (5 q/ 3 in the large-iV approxima- 
tion is found to be 
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Here the spin-orbital index a'^a and we have defined 
several self -energies and functions: The average occu- 
pation number of local electrons at the impurity site is 
defined as 

a + v k a >{c{f a >) 
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and the three self -energies are 
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and 

S3,qq' = \ V ka ' 
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where we have defined two functions: -Di. QQ ' {k,Lu) = 
uj + e k -e a -e a > - U and D 2 , aa ' {k,uj) = u) - e k - e a + e a ,. 

The final expression of Eq. 10 is nontrivial and re- 
quires a delicate derivation. We can simply summa- 
rize our main results: 1) High order Green's functions 
have been naturally included in our derivation, which 
enables us to truncate arbitrary order to close the cou- 
pling integral equations; 2) A higher order decoupling 
scheme we have used in our derivation makes our for- 
mula distinguish from the previous conventional N = 2 
results; 3) Eq. Q is not an extension to APL Ja, |6(], but 
a new large- decoupling scheme. It is due to the fact 
all of the terms N — 1, which are implicitly included in 
Eq. l[3]l and explicitly expressed in Eqs. (4) and (6) , are 
absent if one takes APL's decoupling with large- N EOM. 

We shall illustrate the new features and ingredients 
of Eq. 10: 1) It includes a set of the N>2 closed self- 
consistent integral equations, which can be numerically 
solved; 2) The effective occupancy is frequency depen- 
dent, 3) The higher order self-energy contains the in- 
termediate off-diagonal states in momentum space (e.g., 
(c k c q )) and charge fluctuations (e.g., (f a c q )); 4) It repro- 
duces the results from the standard TV = 2 EOM in the 
infinite-?/ limit H. 

Before we proceed to carry out numerical calcula- 
tions, we shall also point out that the expectation value 
of (c q c k ) and (f a iC k ), which have been discarded in the 
EOM with APL decoupling scheme prove to be very im- 
portant at low temperatures since they diverge logarith- 
mically at the Fermi level as the temperature approaches 
to zero. Their values should be self-consistently evalu- 
ated through the spectral densities between the conduc- 
tion electron and the impurity Green's functions. 
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where /f.d(w) = l/[exp(u;/fc_BT) + 1] is the Fermi- 
Dirac distribution function, and the Green's function 



<c fcff |c^ ff > is: 
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and similarly 

(flc k ) = -- I ' f(w)Jm<c k \fl»dw , (10) 
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Figure 1: (Color) (top): Spectral density DOS calculated via the 
EOM method for particle-hole symmetric Anderson impurity 
model for different orbit-spin degeneracy N at zero tempera- 
ture with Coulomb interaction (7 = 4, solid lines correspond 
to our large N > 2 formulation and dash line is the previous 
TV = 2 formulation result; (bottom) Results of spectral density 
DOS via large TV limit EOM for N — 2 particle-hole symmetric 
Anderson model with different Coulomb interaction U 



Solving the coupled Eq. not only yields the correct 
Kondo resonance at low temperatures but also allows 
us to explicitly include the logarithmic divergence in 
general. When the spin-dependent effect is taken into 
account, the significance of our new EOM approach is 
that it does not rely on the additional renormalization 
introduced in the previous EOM technique fioh . The 
purpose of the additional renormalization is to account 
for the spin dependent level splitting and broadening 
1 23J, 12411 . The lack of rigorous justification for the exis- 
tence of the additional renormalization has cast a doubt 
for the effectiveness of the EOM approach for the non- 
equilibrium Kondo problem. In our improved EOM for- 
mula, we find that the correct Kondo resonance can be 
derived without introducing the additional renormal- 
ization. Comparing with previous calculations, we have 
properly evaluated terms such as (clck), (fl>Ck) through 
Eq. ^ [9]. These terms make crucial contributions to the 
Kondo resonance peak at very low temperatures. Ne- 
glecting these terms will lead to severe errors, which has 
to be recovered by adding ad hoc an additional renormal- 
ization. 



In order to demonstrate the power of our large- N 
EOM approach, we shall apply our new formulation to 



consider the Kondo impurity problem in the particle- 
hole symmetric case, where the standard N = 2 EOM 
fails to show the Kondo resonance peaks. In this case, 
the impurity levels for Vj. Q = are symmetrically placed 
about the Fermi level at e = ep — U/2 and e = ep + U/2. 
If the conduction bands are symmetric about the Fermi 
level and half-filled, the model has complete particle- 
hole symmetry and the average occupation of each spin- 
orbital (n a ) = 1. This symmetric model can display the 
full range of behavior from non-magnetic for ksT,U < 
Aq, to magnetic and Kondo behavior for U^$>Aq, where 
Ao = — Im[So. a (o; + «0 + )]. We shall examine this well- 
studied case by numerically solving Eq. lO and we 
choose the following parameters for our numerical cal- 
culation: The energy of the half-width of the impurity 
resonance in a nonmagnetic metal, Ao is taken 0.1 in the 
unit of conduction band half-width D unless specified 
otherwise. In this special case, the Coulomb interaction 
energy U has usually been taken as a parameter. 

The first illuminating example is the more familiar 
equilibrium Kondo problem where neither impurity 
state nor hybridization is spin-dependence and spin- 
orbital degeneracy N. As shown on the top panel of 
Figure[TJ when one uses the formulation from the stan- 
dard EOM for N = 2 0j, there is no Kondo resonance 
peaks (dashed line). This result numerically confirms 
the analysis made by Kashcheyevs et al. \Vy\ that the 
standard EOM technique, as a severe drawback, can- 
not produce the Kondo resonance at the particle-hole 
symmetric point. Interestingly, within our large- N EOM 
technique, the Kondo resonance peak at the Fermi en- 
ergy is indeed obtained even for N = 2. With increase 
of the spin-orbital degeneracy N, the spectral weight is 
transferred more significantly from the virtual bound 
states toward the coherence region around the Fermi 
energy. This behavior is similar to the case, as demon- 
strated below, when the localized level is moved toward 
the Fermi energy within the energy region of bare reso- 
nance width Ao — mixed valence region. The bottom 
panel of the Fig. [T] shows the impurity spectral density 
(P/T = Pfl) DOS for different Coulomb energy U at zero 
temperature by using our large- N EOM formula. At 
small values of U, the DOS exhibits the two broad peaks 
at e a = —U/2 and e a = U/2, respectively, and a sharp 
Kondo resonance peak at the Fermi level. With the in- 
creased U, the Kondo resonance peak disappears grad- 
ually, evolving into a localized state. These well-known 
results agree with many various approaches, e.g;, the 
scaling analysis jH, the NRG method and NCA j8fl. 

We now consider another very interesting application 
of our large- N limit EOM technique to the Kondo impu- 
rity problem in the presence of external magnetic field. 
We show, in the top panel of Fig. 12 the spectral den- 
sity DOS for a finite U particle-hole symmetric model 
at different exchange fields. A splitting of the Kondo 
resonances peaks for the spin-up and down electrons is 
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Figure 2: (Color) Spectral density DOS calculated via the large 
N limit EOM method N = 2 for a finite U patrticle-hole sym- 
metric Anderson impurity model at different exchange fields 
(top) as well as at various temperatures(bottom). The inset dis- 
plays the zoom-in view of the Kondo resonance near the Fermi 
energy. 



obtained at energy e = ±2H ex . There is no spurious 
peak but only a small bump showing up near the Fermi 
energy. This result is in reasonable agreement with that 
from the NRG calculations l25l l26ll , except for a small 
bump obtained here. The temperature dependent effect 
has been shown in the bottom panel of Fig. [2] . We find 
that the width and height of Kondo resonance peaks 
dramatically changed with increasing temperature. The 
inset displays the zoom-in view of the Kondo resonance 
near the Fermi energy. By further increasing the temper- 
ature, the Kondo resonance peaks disappear at a charac- 
teristic temperature (about 0.3 for the given parameter 
values). We should also mention but not show that the 
new large-TV EOM technique can describe equally well 
the Kondo physics for an asymmetric Anderson impu- 
rity model. 

In conclusion, we have developed a large-A^ EOM 
approach to the Kondo impurity problem for arbitrary 
Coulomb interaction U at finite temperatures. Numer- 
ical results are carried out for symmetric Anderson im- 
purity model with finite U . We show that the Kondo 
resonance peak, which escapes from the standard EOM 



approach for the particle-hole symmetric point, can be 
restored in the new technique. Furthermore, we have 
also shown that the new technique describes reasonably 
well the field dependence of the Kondo effect. Both suc- 
cesses establish the power of this new EOM technique. 



We here should give special thanks to Shufeng Zhang 
and R. C Albers for useful discussions. This work was 
supported by the Robert Welch Foundation No. E-1146 
at the University of Houston (Y.Q. and C.S.T.), by U.S. 
DOE under Contract No. DE-AC52-06NA25396 and un- 
der Grant Nos. LDRD-DR X9GT & X9HH (J.X.Z.). 



* Electronic address: yqi@mail.uh.edu 
t Electronic 



jxzhu@lanl.gov 



[4] 



address: 

URL: |http : / /theory . lanl . gov| 
* Electronic address: csting@mail.uh.edu 
[1] P.W. Anderson, Phys. Rev. 124, 41 (1981). 
[2] A. C. Hewson, The Kondo Problem to Heavy 

Fenwzons(Cambridge University Press, Cambridge, 

U.K., 1993). 

[3] H.R. Krishna-murthy J. W. Wilkins, and K. G. Wilson, 
Phys. Rev. B 21, 1003 (1980); 21, 1044 (1980). 
A.M. Tsvelick and P. B. Wiegman, Adv. Phys. 32, 453 
(1983). 

[5] J. A. Appelbaum and D. R. Perm, Phys. Rev. 188, 874 

(1969) . 

[6] C. Lacroix, J. Phys. F: Metal Phys. 11, 2389 (1981). 
[7] G. Czycholl, Phys. Rev. B 31, 2867 (1985). 
[8] N. E. Bickers, Rev. Mod. Phys. 59, 845 (1987). 
[9] Y. N. Qi, et. al, Phys. Rev. B 78, 045305 (2008). 
[10] Y. Meir, N. S. Wingreen, and P. A. Lee, Phys. Rev. Lett. 66, 
3048 (1991). 

[11] J.A. Appelbaum and D. R. Perm, Phys. Rev. B 3, 942 (1971). 

[12] A. Theummann, Phys. Rev. 178, 978 (1969). 

[13] H. Mamada and F. Takano, Prog. Theor. Phys. 43, 1458 

(1970) . 

[14] G.S. Poo, Phys.Rev. B 11, 4606 (1975); 11, 4614 (1975). 
[15] V. Kashcheyevs et al, Phys. Rev. B 73, 125338 (2006). 
[16] Th. Pruschke and N. Grewe, Z. Phys. B 74, 439 (1989); K. 

Haule et al, Phys. Rev. B 64, 155111 (2001). 
[17] P. Coleman, Phys. Rev. B 29, 3035 (1984). 
[18] N. Read and D.M. Newns, J. Phys. C 16, L1055 (1983). 
[19] EC. Zhang and T.K. Lee, Phys. Rev. B 28, 33 (1983); ibid. 

30, 1556 (1984). 

[20] J. W. Rasul and A.C. Hewson, J. Phys. C 16, L933 (1983). 
[21] J. W. Rasul, N. Read and A.C. Hewson, J. Phys. C 16, 
L1079 (1983). 

[22] F. D. M. Haldane, Phys. Rev. Lett. 40, 416 (1978). 
[23] Y. Meir, N. S. Wingreen, P. A. Lee, Phys. Rev. Lett. 70, 2601 
(1993). 

[24] J. Martinek et al, Phys. Rev. Lett. 91, 127203 (2003). 
[25] T. A. Costi, Phys. Rev. Lett. 85, 1504 (2006). 
[26] A. C. Hewson, J. Bauer, and W. Koller, Phys. Rev. B 73, 
045117 (2006). 



